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Abstract. We report results of a new analysis for the low-temperature structure of a'-NaV20s from syn- 
chrotron x-ray diffraction experiments. We confirm the existence of two inequivalent ladder structures in 
each vanadium layer. Based on our structural data we perform a bond- valence calculation for the vanadium 
sites in the low temperature state. Due to an asymmetric charge ordering we obtain only two different 
vanadium valences despite the three inequivalent sites. This explains the 51 V-NMR observation of only two 
resonant peaks in the charge ordered phase. By use of a Slater-Koster method to obtain hopping matrix 
elements and cluster calculations we obtain effective vanadium- vanadium hoppings which compare well to 
LDA results. Using these in a cluster calculation we obtain a superexchange of 0.047 eV between electrons 
on neighbouring rungs of the same ladder for the undistorted phase. For the distorted phase we find a 
significant alternation in the shifts of the oxygen atoms along the legs of one of the two ladder types which 
leads to a significant exchange dimerisation 8j « 0.25. 



PACS. 61.66.Fn Inorganic compounds - 61.50.Ks Crystallographic aspects of phase transformations; 
pressure effects - 75.30.E Exchange and superexchange interactions 



1 Introduction 

The layered oxide cv'-NaV2 05 has attracted great interest 
since 1996, when Isobe and Ueda reported a phase transi- 
tion at T = 34K with a spin-Peierls like spin gap forma- 
, tion jp . An earlier determination of the high-temperature 

■ structure 0] originally reported two inequivalent vana- 
dium sites. However, a redetermination of the crystal struc- 
ture of the high-temperature phase found the space group 
rmrnn |,|,§ instead of P2 imn, with only one type of 
vanadium site. Due to stoichiometry this site should have 

■ a valence of +4.5. Later, a 51 V-NMR study by Ohama et 
al H reported the existence of two inequivalent vanadium 
sites below the critical temperature. This is supposed to 
be caused by a charge ordering transition of the mate- 
rial which leads to inequivalent vanadium atoms with dif- 
ferent valences. A determination of the low-temperature 
structure by x-ray scattering then reported j^i/iree in- 
equivalent vanadium sites. Using this data a bond-valence 
calculation then led to three different valences, B approxi- 
mately given by +4, +4.5, and +5, in contrast to the result 
of the 51 V-NMR measurements. 

In this article, we will give the results of a determina- 
tion of the low-temperature structure using a much larger 
data set than used in Q. Applying bond- valence analy- 
sis to our data we obtain only two significantly different 



valences for the vanadium sites, despite the fact that our 
structural data show also three inequivalent sites. 

Due to shifts in the oxygen positions we further find 
an alternation in the exchange coupling constant between 
the vanadium sites on one type of ladder and give an esti- 
mate for the corresponding exchange dimerisation param- 
eter 5 j. 

This article is organised as follows. In the next sec- 
tion we present the experimental details and the results 
of the structure determination. We find three inequiva- 
lent vanadium sites in each layer, building two inequiva- 
lent vanadium-oxygen ladders per layer. In the third sec- 
tion we use the bond-valence method to determine the 
valences of the sodium and the vanadium atoms both in 
the undistorted and the distorted phase and compare the 
results with the observations of NMR. In the fourth sec- 
tion we look at the superexchange coupling along ladder 
direction for the undistorted phase and for the two ladders 
per layer of the distorted phase. To do this we first deter- 
mine approximate vanadium-oxygen and sodium-oxygen 
hopping matrix elements. From these we find the effective 
vanadium-vanadium hopping elements by use of a block 
diagonalization comparing our results with those of a re- 
cent LDA calculation [Qj . We proceed to find the superex- 
change coupling for the different ladder types and deter- 
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Table 1. Average structure parameters obtained by using fundamental reflections at 15K. 



Site x y z Site 

V 0.40212(7) 0.25000(0) 0.39067(22) Na 

Oi 0.25000(0) 0.25000(0) 0.51111(122) O a 

3 0.42595(39) -0.25000(0) 0.50381(93) Cell [[| 



mine its alternation along the ladder direction. In the last 
section we summarize and discuss our results. 



2 Structure details 

The crystal structure of the room temperature phase of 
NaA^Os based on X-ray diffraction intensity data has 
been reported previously || . We have measured the diffrac- 
tion intensities of NaV205 at T = 15 K in the charge or- 
dered phase using the same crystal. The experiment was 
performed using a single crystal of size 0.03 x 0.09 x 0.20 
mm 3 on the diffractometer D3 at the synchrotron X-ray 
source of the HASYLAB at DESY, Hamburg by using 
X-ray wavelength A = 0.559A. The sample was mounted 
inside a Displex cooling device attached to the four-circle 
diffractometer. The superlattice reflections were detected 
corresponding to the lattice parameters a = 2 x ao, b — 
2 x bo and c = 4 x cq where ao, bg and Co are subcell lat- 
tice parameters. To determine the space group, reflections 
were measured using the primitive super cell for small val- 
ues of Bragg angles which showed that the space group is 
Fmm2. 4847 reflections were measured at T= 15K up to 
sm 9/X = 0.94A- 1 which gave rise to 1953 unique reflec- 
tions. We have determined the average structure by using 
only the fundamental reflections and the room tempera- 
ture structure model in space group Pmmn. Least-squares 
refinement of the crystal structure by using SHLEX pro- 
gram led to the atom parameters listed in table [|. The 
residual index R, defined as 

rt Y.\\Fo\-\F c \\ ( 

was R — 0.060. The average structure obtained agrees 
closely with the room temperature structure^ as expected. 
We next refined the charge-ordered distorted low temper- 
ature structure in the space group Fmm2. The number of 
parameters refined were 61. The starting parameters were 
the parameters of the average structure (table |l|) trans- 
formed appropriately to the supercell. The resulting atom 
coordinates have been given in table [|. Note that contrary 
to the superspace group refinement of reference we 
have employed the standard crystallographic refinement 
procedure which we could afford due to the very large 
number of reflection intensities measured up to a large 
value of sin 9/ A = 0.94A -1 . Accordingly our parameters 
were free from the constraints of the refinement using the 
superspace group. Our refinement also led to physically 
meaningful thermal parameters. The residual index was 
R= 0.056. The atom coordinates are listed in table g. To 



x y z 

0.25000(0) -0.25000(0) 0.85362(73) 

0.38565(42) 0.25000(0) 1.05413(151) 

11.3030(1)A 3.61095(3)A 4.7525(1)A 



ensure that the refinement did not lead to a local minimum 
we used the so-called multistart-and-refine procedure in 
which one starts from several initial displaced coordinates. 
These refinements lead to identical coordinates within the 
standard deviations. Recently doubts have been raised 
about the the actual space group of NaV^Os. 23 Na-NMR 
results have suggested a lower-symmetry space group than 
Fmm2 |l(| . We did not attempt to refine the crystal struc- 
ture by lowering the space group symmetry because of the 
very large number of parameters and consequent strong 
correlations among the parameters involved. Also there 
exists no diffraction evidence for a lower symmetry space 
group. Details of the present refinement procedure will be 
published elsewhere pl[ . 

We get the same qualitative results as were found in (?], 
|l2f . Two slightly different types of vanadium layers were 
found. Both layers look similar as shown in figure [I] with 
three inequivalent sites and as described in Q] . We denote 
the two different layers by the subscripts a and b and 
the vanadium sites as shown in figure nL Apart from the 
slightly different magnitude of the shifts the main quali- 
tative difference is, that for layer a the Vi a sites shift in b 
direction towards the V2i a sites while for layer b the Vk, 
sites are found to shift away in b-direction from the V2ib 
sites. 



3 Bond-valence calculation 

The bond- valence method serves to determine the valences 
of atoms in a chemical compound jl4j . The valence of an 
atom in a compound, as determined by this method, is 
the sum of the bond valences for all bonds, in which the 
atom participates. The bond valences of such anion-cation 
bonds are empirically defined as 

Vi = exp [(ro - n) I B] , (2) 

where B is fixed to be 0.37. r, is the distance between 
the two atoms of the anion-cation bond, ro is an empirical 
parameter obtained from fitting the atomic valences for 
many different compunds at room temperature [fl5| . There 
is an alternative definition for the bond valence, namely 




(3) 



where both ro and N are empirical parameters. Both def- 
initions give similar results; we will adopt definition u2) 
throughout this article, as this method was also used in ||. 

It should be noted that the bond valence model does 
not imply a certain physical nature of the bond. It works 
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Table 2. Low temperature phase structure at T = 15 K and shifts from high temperature structure as given in |3j]. Size of 
unit cell taken from Shifts from the positions in [^) have been calculated after shifting the low temperature unit cell by 
the vector (-1/8, 1/8, 0.87959/4). The shift in z direction is chosen, such that the positional change of V2 and V4 is nearly zero 
and symmetrical. Thus, the vanadium sites of the spin dimerised ladders do not move in z direction below Tc- 
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Fig. 1. Vanadium layer of the low-temperature structure. Note the difference in the shifts of Via and Vi6- Figure not to scale. 
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Fig. 2. Oxygen (white) and vanadium atoms (grey) around 
vanadium atom (black) on same ladder. Arrows indicate shift 
of oxygen atoms along the leg in the Vi-ladders. Figure not to 
scale. 



both with covalent and ionic type bonds Jig ]. Inaccura- 
cies generally come from two sources: First, the errors of 
the structure determination and the lattice constant de- 
termination. Second, the inaccuracy of the parameter r§ 
in (0). This parameter is determined from fitting it to a 
large number of structures containing a given bond type 
and therefore carries an error with it. We give the total 
standard deviation derived from these two sources of er- 
rors in our results for the valences. In our calculations 
for the vanadium site valences the inaccuracies from the 
structure determination and the lattice constant play the 
leading role. 

Therefore, in a first step we use the method to deter- 
mine the valences of the atoms in the undistorted phase. 
This enables us to adjust parameters to account for the 
peculiarities of the material as well as to estimate the ac- 
curacy of the method. In a second step we determine the 
valences of the vanadium and sodium atoms in the dis- 
torted phase. This will give us some insight into likely 
charge order structures. 



3.1 Atomic valences in the undistorted phase 

We have obtained the bond lengths for the vanadium- 
oxygen bonds and the sodium-oxygen bonds from the struc- 
tural data in |1 and give them in tables |^ and ^. 

Using the lengths of the vanadium-oxygen bonds we 
determine the valences of the vanadium and the sodium 
atoms in the undistorted phase by Using room tem- 
perature lattice parameters and the empirical tq given 
in jTa we obtain two different results, depending on which 



Table 3. Vanadium-oxygen bond lengths given in An gstrom. 
Vo denotes the vanadium site of the undistorted phase. For 
the nomenclature of the oxygen sites, see figure bl The first 
row gives the distances for room temperature lattice parame- 
ters as given in the second row gives distances for lattice 
parameters at 15K — see table ^. 





Oi 


o 2 


o 3 


o 4 
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Vo 
Vo 


1.916(3) 
1.914(3) 
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1.8216(7) 


1.916(3) 
1.914(3) 


1.986(3) 
1.983(3) 


1.616(3) 
1.600(3) 



Table 4. Sodium-oxygen bond lengths given in Angstr0m for 
room temperature lattice constants (first row) and lattice con- 
stants at 15K (second row). NaO denotes the sodium site of the 
undistorted phase; for the nomenclature of the oxygen sites see 
figure []. Due to symmetry there are always pairs of bonds with 
the same length, thus, the lengths of only four bonds are given. 
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NaO 
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value we take for ro. Using the value for V 4+ -0 bonds 
i'4 + = 1.784(27)A gives a valence of 4.45; using the value 
for V 5+ -0 bonds r5 + = 1.803(31)A leads to a valence 
of 4.69. We know from stoichiometry that the vanadium 
atom should have a valence of 4.5. Just using the aver- 
age of the two ro given above is not a good approxima- 
tion, since the characteristic bond lengths do not vary 
linearly with the bond order. Therefore, we use a fitted 
r 4 . 5+ = 1.788(5)A to give the correct valence 4.5 for the 
undistorted phase with room temperature lattice param- 
eters. 

The size of the unit cell decreases significantly from 
300K to 15K. We therefore cannot simply use the r<i + , 
r4.5 + and r§ + from above, since they result from room tem- 
perature structure determinations. Using the high temper- 
ature r4.5 + with the low temperature lattice parameters 
for determination of the valence of the vanadium site in 
the undistorted phase would give a valence of 4.59. We 
therefore have to scale down all these characteristic bond 
lengths. For this scale we use the factor (V15K /Vzook) 1 ^ = 
0.9965 where the V are the volumes of a cell containing 
four vanadium sites. We then obtain a valence of 4.51 for 
the vanadium atom in the undistorted phase using low 
temperature lattice constants and the reduced r4.5 + . 

We can follow through a similar exercise for the sodium 
atom. We find that the characteristic bond length r^ a = 
1.803A given in [|l5| gives a valence of 1.12 for the sodium 
atom in the undistorted phase with room temperature lat- 
tice parameters, which is clearly deviating strongly from 
the expected value 1. It is known, however, that there 
is a comparatively strong physical variation in the envi- 
ronment of alkali atoms, therefore the characteristic bond 
lengths vary strongly, the standard deviation being about 
0.08A |lq| . It is therefore justified to use a characteris- 
tic bond length r^ a = 1.762(5)A obtained by fitting the 
resulting valence to 1. This helps capturing the special 
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Fig. 3. Sodium (black), oxygen(white) and vanadium (grey) 
atoms on two ladders in two neighbouring layers. Figure not 
to scale. 

physical features of the environment surrounding the al- 
kali atom in this material. Scaling down this character- 
istic r^a to 1.756(5) for lattice constants at 15 K as in 
the case of vanadium and using these lattice constants we 
arrive at a valence of 1.00(1) for the sodium atom in the 
undistorted phase using low temperature lattice param- 
eters. Since r/y a is a fitted value, the error comes only 
from the standard deviation of the structure determina- 
tion and is much smaller than the standard deviation of 
TMa from . All the parameters used are listed in table 0. 
Henceforth we will always use the lattice parameters and 
characteristic bond lengths for temperature 15K as given 
there. 

To check the quality of the bond valence method we 
now determine the valence of the three oxygen sites of the 
undistorted phase. Since the bond orders are already de- 
termined, the total valence of the five oxygen atoms in a 
unit cell is —10.02(6), corresponding to the total valence of 



Table 5. Parameters used for the determination of atomic 
valence. Lattice constants from j|[l3|], for other parameters 
see text. 
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the cations. However the valence of the different oxygens 
is not the nominal value —2 for each as might be expected. 
Instead we obtain a valence of —1.87(2) for the apical oxy- 
gen, —2.08(2) for the oxygen on the leg of the V-O-ladder, 
and —2.12(2) for the oxygen on a rung. Qualitatively this 
may be explained with the fact that the oxygen atoms are 
easily polarizable. 

3.2 Atomic valences in the distorted phase 

We will now determine the valence of the vanadium and 
the sodium atoms in the distorted phase at 15K. For this 
we will use parameters as given in table [s]. 

The bond lengths of the vanadium-oxygen bonds and 
of the sodium-oxygen bonds as well as the valences result- 
ing are listed in tables ^ and 0. Due to the fact that the 
exact valence of the vanadium atoms is not known a pri- 
ori and that the characteristic bond length ro depends on 
this valence, one has to use a self consistency procedure 
to derive the valences for the vanadium atoms as given in 
the last column of table ^. 

Regarding the sodium atom valence one observes, that 
the changes of valence are smaller than the standard de- 
viation, which is mostly due to the standard deviation 
of rjva- Furthermore the average valence does not change 
significantly, it is equal to 1.01(3). 

Regarding vanadium, the first observation is that the 
valences of similar vanadium atoms in the two layers are 
not the same within one standard deviation. We will there- 
fore look at these layers separately. 

We first look at the valences for layer a. Although the 
Vi a and the N^ia sites are crystallographically inequiva- 
lent they have approximately the same valence 4.6. The 
V2ia site has a lower valence of approximately 4.1. Note 
that the shift from the average valence 4.5 is significantly 
different for V2i a and V22a- 

For layer b the results are somewhat different. Here 
the Vib site experiences only a small change of valence. 
Despite this the absolute size of the valence shifts from 
the original average valence 4.5, which are 0.43(4) for V216 
and 0.19(4) for V22&, are significantly different from each 
other. 

To determine the quality of the results, one can cal- 
culate the stoichiomctrically averaged valence of a vana- 
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Table 6. Vanadium-oxygen bond lengths given in Angstr0m and valences in the distorted phase at 15K. Nomenclature for 
atoms as given in figures g and |l| For the valences The first three columns show the valence depending on the characteristic 
bond length ro used, the last row gives values which are self-consistent to first order. 
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4.57(4) 


4.76(4) 


4.61(4) 


Vlf, 


1.914(3) 


1.821(4) 


1.909(3) 


1.992(7) 


1.608(7) 


4.43(4) 


4.48(4) 


4.66(4) 


4.48(4) 


V 2 1 6 


1.938(3) 


1.900(7) 


1.938(3) 


1.981(7) 


1.635(6) 


4.07(4) 


4.11(4) 


4.28(4) 


4.07(4) 


V226 


1.892(3) 


1.752(7) 


1.892(3) 


2.004(7) 


1.632(6) 


4.57(4) 


4.62(4) 


4.81(4) 


4.69(4) 



Table 7. Sodium-oxygen bond lengths given in Angstr0m and 
valence in the distorted phase at 15K. Nomenclature for oxygen 
atoms as shown in figure []. Due to symmetry there are pairs 
of bonds of equal length, only one of each pair is given. 





Or 


Ol 


Oai 


OA2 


Valence 


Nal 


2.405 


2.570 


2.513 


2.550 


1.06(8) 


Na2 


2.418 


2.583 


2.541 


2.515 


1.04(8) 


Na3 


2.444 


2.569 


2.576 


2.576 


0.97(8) 


Na4 


2.443 


2.568 


2.567 


2.567 


0.98(8) 


Na5 


2.423 


2.630 


2.572 


2.572 


0.96(8) 


Na6 


2.420 


2.621 


2.548 


2.548 


1.00(8) 



dium site for the low-temperature structure. Due to charge 
conservation one would ideally expect an average valence 
of 4.5. Using the self-consistently interpolated values one 
obtains 4.48(2) for layer a and 4.43(2) for layer b. This 
indicates that the results for layer a should be the more 
reliable. The difference between the results for the differ- 
ent layers can well be caused by the refinement method 
of the structure determination, which might have found 
a less reliable local minimum for layer b than for layer a. 
From these bond valence calculations for layer a we there- 
fore conclude that there are only two significantly differ- 
ent vanadium valences in the distorted phase, contrary 
to the result of van Smaalen and Liidecke, who obtained 
three clearly different valences Q and a nearly symmetri- 
cal charge shift between the V21 and the V22 sites. How- 
ever, in our analysis we obtain a significantly asymmetrical 
charge shift between the V21 sites on one hand and the 
Vi and V22 sites on the other hand for both layers. 

Such a charge transfer between ladders does not have 
to cause the material to become a conductor. The insulat- 
ing behaviour of the high temperature phase has been ex- 
plained with correlation effects, due to which one can map 
the quarter filled ladder to a half-filled Hubbard chain |t7| , 
resulting in an effective charge transfer gap. Since this ar- 
gument works only at quarter-filling, one might expect the 
material to become a conductor if a charge transfer takes 
place between ladders as described above. While this is 
true if one considers isolated ladders, the situation be- 
comes less simple if inter-ladder Coulomb repulsion is in- 
cluded in the model, due to the geometry of the lattice. 
Consider a hexagon of Vi and V22 sites and assume that 



the V21 sites carry one electron each. In this case one has 
two electrons for each six other sites, close to what is found 
from bond valence above. As a consequence, one still has 
a charge transfer gap, because due to the geometry of the 
lattice it is not possible to put three electrons on such 
a hexagon without having Coulomb repulsion. Therefore 
depending on the ratio of the parameters involved, the 
system can still be a charge transfer insulator. 



3.3 Comparison with NMR 

We can compare our result to the results of the 51 V-NMR 
measurements ||. For such a measurement one expects 
that the resonant peak for the high-temperature vanadium 
site with valence +4.5 splits into multiple peaks below the 
transition: one peak for each different valence of vanadium 
sites is expected. The intensity of a peak for a given va- 
lence should be approximately proportional to the number 
of vanadium sites having such valence. 

The interaction between the 3d electrons and the nu- 
clei of the vanadium which causes the shift of the NMR 
resonance frequency, is, disregarding crystal field effects 
due to the distortion, roughly proportional to the den- 
sity of the d^y-electrons on a shell of a given vanadium 
ion. Therefore one expects that the splitting of the reso- 
nant NMR-peak at the transition temperature reflects the 
charge disproportionation approximately linearly, at least 
close to the critical temperature. 

For a low-temperature structure with three different 
valences 4.5 — S for the V21 sites, 4.5 for the Vi sites, and 
4.5 + 5 for the V22 sites, as suggested by one expects a 
splitting of the single high-temperature resonant peak into 
three different peaks with an intensity ratio of approxi- 
mately 1:2:1, the corresponding shifts of the right and 
left peak should be roughly symmetrical around the cen- 
ter peak and the center peak in continuation of the single 
high temperature resonance peak. For a low-temperature 
structure with only two different valences 4.5 — 3<5 for the 
V21 sites, and 4.5 + 5 for the V22 and the Vi sites as sug- 
gested by our analysis above, one expects a splitting of 
the high-temperature resonant peak into two peaks at the 
transition temperature with an intensity ratio of approx- 
imately 1 : 3 and where the shifts of the new peaks from 
the old peak obey a 3 : 1 ratio. 
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Looking at the experimental results as given in figure 

3 of ||, one finds, that two resonant peaks have been 
found below 33. 4K. The shift of these two peaks relative to 
the high-temperature peak is clearly asymmetric, and the 
ratio of the distances is roughly 2 : 1 in favour of the left 
peak which has been identified to be a V 4+ resonant peak 
by the authors. The ratio of the intensities, as determined 
by the area below the peaks, is roughly 1 : 3 from 30K 
to 33. 3K. This changes at 20K, but it is not mentioned 
in |(| whether intensity correction has been performed for 
the data, so as far as intensity ratios are concerned we use 
only the data around Tc- 

From 33. 4K to 33. 6K three peaks are found. One of 
them is clearly the continuation of the high-temperature 
peak for valence 4.5. This behaviour might be intrinsic 
or due to a distribution of the Tc m the powder sample 
used for the measurement. In the former case (intrinsic 
behaviour) one would identify the center peak with the 
Vi sites and the left and right peak with the V21 and the 
V 2 2 sites, respectively. Since the charge of the Vi sites 
would then remain constant in this temperature interval, 
the charge shift of the V21 and the V22 sites would be of 
the same amount and different sign. This implies a ratio 
of roughly 1 : 1 for the distances of the two side peaks 
from the center peak, i.e., a different ratio than that mea- 
sured below 33. 4K. For the latter case (distribution of Tc) 
one would not expect a big change of the distance ratio 
at 33. 4K. Since one finds that the ratio of distances is 
still 2 : 1 from 33. 4K to 33. 6K, we believe that the ob- 
servation of three peaks is due to a distribution of the 
Tc rather than intrinsic. However, below 33. 4K the 51 V- 
NMR measurements seem to be in much better agree- 
ment with the results of our analysis, yielding only two 
different valences with an asymmetrical charge dispropor- 
tionation, than with the previous analysis of van Smaalen 
and Liidecke yielding three different valences and sym- 
metrical charge disproportionation Our results also 
explain, that the low-valence V21 sites can have a valence 
of nearly 4 as determined in Q from the hyperfine cou- 
pling strengths, although the charge disproportionation Sq 
in the system is estimated to be small: 8q <C 0.5 |T8|.|r9| . 

4 Exchange coupling constant alternation 

From the structural data we observe a slight alternation 
in the distances in b-direction between Vi sites for the Vi 
ladder in both layers. The lattice dimerisation parameter 
Sy for the shift of the Vi sites in b direction is very small 
with 0.04% for layer a and 0.05% for layer b compared 
to 0.56% for CuGeOs |20]| . However, even a small alter- 
nation implies an alternation in the exchange coupling in 
b-direction. 

Under the assumption that the observed spin-gap re- 
sults from an alternation of the exchange coupling in all 
ladders, the exchange dimerization parameter Sj = (J+ — 
J_)/(J+ + J ) is estimated to lie between 3.4% to 10% for 
a'-NaV 2 5 [£||§|||| and 3% to 17% for CuGe0 3 fcj, 
|25j . In both cases the higher values result from mean- field 
estimates, which might be insufficient for a treatment of 



these materials as argued for a'-NaV205 in |2]j| . However, 
the assumption that there is an exchange coupling alter- 
nation in all ladders of the distorted phase is flawed from 
the point of structure determination, since we observe only 
for the Vi-ladders any alternation in b-dircction, but not 
for the V2i-V22-ladders. 

Assuming that an alternation of the exchange coupling 
in b-direction takes place only in every second ladder, a 
recent theoretical study estimated the dimerisation pa- 
rameter S necessary to produce the observed size of the 
spin gap as 5 » 0.38 . . . 0.54. Q This would be even 
larger than the dimerisation observed in CuGeOs. Thus, 
it seems unlikely, that the shifts of the vanadium sites and 
the resulting dimerisation in a'-NaV2 0s which is an order 
of magnitude smaller than that for CuGeOs should give 
an exchange dimerisation of the same order of magnitude. 

However, the small dimerisation of the vanadium sites 
is not the only alternation along the legs of the Vi ladders. 
One finds a much stronger alternation in the distances to 
the nearest oxygen atoms which are on the same leg as 
the Vi atom in question. This alternation results from 
the shift of these oxygen atoms in a- and c-direction. For 
the V21-V22 ladder no effective alternation in b-direction 
is found, neither from shifts of the oxygen sites nor from 
shifts of the vanadium sites. 

It is possible to determine the exchange dimerisation 
parameter from the structural data. Let us first give a very 
simple estimate. For the superexchange interaction be- 
tween rungs we have approximately J cx t vv , tyy oc t vo , 
and, according to [^7), tyo d v ^f. tyo is a hopping ma- 
trix element between the vanadium d^-orbital and the p- 
orbital of the neighbouring oxygen on the same rung, tyy 
is therefore an effective hopping matrix element between 
neigbouring vanadium sites of on the same leg of a lad- 
der. This implies, that for small distortions which change 
the distance dyo by ±Adyo we obtain Sj = AJ / J w 
±\AAdvo/dvo- For the vanadium-oxygen bond lengths 
between atoms on the same leg of the Vi ladder given in 
table [| one finds 5j w 2%. This is one order of magni- 
tude larger compared to what one would have obtained by 
using the shifts of the vanadium sites only. 

We will now first examine the hopping matrix ele- 
ments for the undistorted phase, calculating them with the 
method presented in psfl using the parameters from p7[ . 
From this we will gain insight in the effective tyy hopping 
matrix elements and the effective superexchange between 
rungs on the same ladder. In a second step we extend this 
method onto the distorted phase to find the approximate 
size of the superexchange dimerisation on the Vi ladders. 

4.1 Hopping and superexchange in the undistorted 
phase 

We determine the matrix elements for the electron hop- 
ping by a Slater-Koster type of treatment presented in [ p7| 
and fl28fl . We will use hopping matrix elements along anion- 
cation bonds and between vanadium orbitals only. For 
the vanadium sites, we consider only the lowest lying 3d- 
orbital, which is the d xy orbital according to ||. For the 
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Table 8. Distances (in Angstr0m) between vanadium and oxygen atoms and hopping matrix elements (in eV) between vanadium 
d xy and oxygen p-orbitals for the undistorted phase. The absolute values of the hoppings and the distances are given, the sign 
is determined by the orientation of the vector between the vanadium and the oxygen site. Denomination of oxygen atoms as in 
figure 0. O5' denotes the apex oxygen of the VO5 pyramid lying above or below the central vanadium site. 





Z\x 




Az 


dvo 






^Pzd xy 


Oi, 3 


0.2825 


1.8055 


-0.5695 


1.914 


0.8549 


0.6303 


0.2470 


o 2 


1.7193 





0.6018 


1.822 





1.1632 





4 


1.9304 





0.4545 


1.983 





0.8909 





0, 


0.1855 





1.5896 


1.600 





0.2247 





5 ' 


0.1855 





3.1629 


3.168 





0.0104 






Table 9. Distances (in Angstr0m) between sodium and oxygen atoms and hopping matrix elements (in eV) between sodium 
3s and oxygen p-orbitals for the undistorted phase. The absolute values of the hoppings and the distances are given, the sign is 
determined by the orientation of the vector between the vanadium and the oxygen site. Denomination of oxygen atoms as in 
figure |. 





Ax 




Az 


djVaO 




tspy 


^>sp z 


Oai, Oa2 
Ol 

Or 


1.5338 
2.0018 



1.8055 


1.8055 


0.9434 
1.6501 
1.6178 


2.550 
2.594 
2.424 


1.2970 
1.6076 



1.5267 


1.7768 


0.7977 
1.3252 
1.5921 



The coordinate differences, distances and resulting hop- 
ping elements between the sites are given in tables || and |^. 
For the direct d xy -d xy hopping matrix element one obtains 
0.2231 eV in the undistorted phase. 

Using these matrix elements, we will now proceed to 
determine effective vanadium-vanadium hopping in the 
following way. Consider the single particle Hamiltonian 
Hsp of the problem. In matrix notation it has the form 

(Tyv Tvo TvNa \ 
Toy Too T OA f a . (6) 
T^NaV TjVaO TAT a jv a J 

The on-site energies ey, eo, an d fjv a on the diagonal of 
T V v, Too, and T NaNa are given in |^7| as ey = —12.55 
eV, eo = —14.13 eV, ejv a = —5.13 eV. One observes, that 
the hopping from one layer to another via the interme- 
diate apex oxygen is very small. However, this is not so 
clear for hopping via the path V-Oi-Na-CvV. Therefore 
to find effective vanadium-vanadium hopping matrix ele- 
ments we block-diagonalize H given above with successive 
Jacobi transformations, projecting out the sub-matrices 
Tvo, Tatqo and their conjugates to separate the Hilbert 
space of the vanadium orbitals from that of the rest. We 
then proceed to use these effective vanadium-vanadium 
hoppings as parameters of an extended Hubbard model 
describing the system. 

The block diagonalization is done for a cluster of 4 
layers by 4 rungs by 4 ladders. The resulting values for 
the effective hoppings can be found in table HO and we can 
compare them with the results from an LDA tight-binding 
fit presented in 0|. The results are partly similar: t# is 
larger than both ti, and to, although all these values are 
smaller than those from LDA. tz, and ta are of similar size. 
t c , the effective hopping between layers, is much smaller 
than all other effective hoppings in both cases. 



oxygen sites we use the 2p-orbitals, assuming them to be 
of equal energy, while for the sodium sites we use the 3s- 
orbital. In this case, the hopping between vanadium and 
oxygen sites is p8|,|27|] 



t 



m (V?,l 2 V pda + (1 - 2l 2 ) V pd7r ) 
dxy = I (VHm 2 V pdry + (1 - 2m 2 ) V pd ^ (4) 



p < 



with I, m, n being the direction cosines of the vector from 
the p state to the d state and the parameters 

V pda = -21.808 eVA 7/2 dy 7 ^ 2 , 
V pd7T = 10.054 eVA 7/2 dy 7 ^ 2 . 

The matrix elements for the Na-0 hopping are determined 
from 



V spa = 14.0208 eVk 2 d N 2 a0 . 

Finally we may also consider hopping between nearest 
neighbour vanadium sites on different rungs. In this case 
the matrix elements are given by 

td xy d xy = H 2 m 2 V dda + (I 2 + m 2 - 4/ 2 m 2 ) V dd7I (5) 

with parameters 

V dda = -116.18 eVAV 5 



V dd7I = 62.754 eVAV 5 



vo> 
vb- 
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Table 10. Effective vanadium-vanadium hoppings for the undistorted phase given in eV. Denomination of hoppings as given 
in figure 111. For methods see text. Superexchange J e // has been calculated with U = 4 eV, and Vl = 2tft as given (see text). 



Method 


tii 


t L 


to 


tiL 


t c 


Vl 




LDA@ 


0.38 


0.085 


0.085 






0.76 


0.070 


Ab initio |g9j 


0.5382 


0.1246 




-0.0442 




1.076 


0.038 


All hoppings 


0.172 


0.049 


0.062 


-0.110 


0.009 


0.344 


0.047 


Without Na-O 


0.305 


0.100 


0.016 


-0.111 


< 0.001 


0.610 


0.044 



We can also compare our values with the results from 
an ab-initio calculation usin g th e DDCI2 quantum chem- 
ical method as presented in |2{j]. They obtain a t# which 
is larger than ours or that from LDA calculation. They 
do not obtain a value for to because Na orbitals were not 
included in their calculation. Therefore one should com- 
pare their value for tL with our result for trj+tz,. These 
values are similar. They also obtain a hopping matrix el- 
ement tiL which is smaller than ours by about a factor 
of 2.5. These differences are probably due to the fact that 
they include the on-site repulsion energy U p for the oxy- 
gen orbitals. This reduces the energy difference between 
occupation of vanadium and oxygen orbitals and thus in- 
creases the effective hopping between vanadium orbitals. 

A strength of our method is that we are able to analyze 
the originof particular matrix elements. For example, re- 
garding the similar size of tp and tz, we find, that it is due 
to the inclusion of the Na-0 hopping elements. Without 
these we get results as given in the last row of table [To], 
where tjj is lower than tz by nearly two orders of magni- 
tude. Similarly, we notice from table |^ that in our model 
we have more vanadium-oxygen matrix elements than as- 
sumed in the tight-binding fit in As we will see later, 
this is very important for the determination of the effect 
of small distortions. There one can no longer use an effec- 
tive one particle Hamilton operator where only a certain 
subset of all hoppings is included. 

We will next determine the effective superexchange 
J e f / between electrons on neighbouring rungs of the same 
ladder. We do this by considering the Hamiltonian 

+ Vij (n lT + mi) (n n + riji) , (7) 



Regarding U, it has been estimated in || to be 4.0 eV 
and we will use that value. The intersite Coulomb inter- 
action can be estimated in the following way. As shown 
in |3(| the Hamiltonian for the charge degrees of freedom 
can be mapped to an antiferromagnetic fsing model in a 
transverse field, ft is well known [pl| that in such a model 
there is a quantum critical point at 2tn = Vl and it has 
been argued that the material is close to this point |^2[ . 
Using e?z = 3.611 A and <1r = 3.439 A we can estimate 
Vr by assuming that the Coulomb interactions scale with 
due to screening effects. 

The results of our calculation are found in table 10 
The estimate from experiment is 0.045 ... 0.048 eV [33 
|l| , which agrees with the results of our calculation. Using 
the LDA hopping parameters |j) we find a value for J e f f 
from the cluster calculation which is too high. Using the 
parameters from the ab initio calculation |29[ ] we obtain 
a J e ff which is slighly too low, but better than LDA. 
Thus, the set of hopping parameters obtained by the above 
treatment might be preferable for model calculations on 
the ladders than using the LDA results, especially since 
the qualitative features regarding and tz, are conserved. 



where the first term contains the effective hopping be- 
tween vanadium d xy orbitals as obtained either from LDA 
or our method. The second term contains the on-site Coulo 
interaction U, the last term contains the intersite Coulumb 
interaction Vij between sites i and j. To determine the 
effective superexchange we solve this Hamiltonian numer- 
ically on a cluster of two neigburing rungs on the same 
ladder with two d-electrons and set J e f f to be the singlet- 
triplet gap energy on this cluster. We disregard interladder 
hoppings and we will also disregard the small interlayer 
hoppings t c . Now we need values for the onsite Coulomb 
interaction U and the intersite Coulomb interactions. 



4.2 Hopping and superexchange in the distorted phase 

We now repeat the above calculations for the distorted 
phase: We obtain vanadium- oxygen, oxygen-sodium and 
vanadium-vanadium hoppings by using the Slater-Koster 
method. We then project out the vanadium-oxygen and 
the vanadium-sodium hoppings in the single-particle Hamil- 
tonian to obtain the effective vanadium-vanadium hop- 
pings. We calculate the effective superexchange J e // for 
the Vi-ladders as we did for the undistorted phase. How- 
ever, we now have two inequivalent two-rung-clusters on 
this ladder due to distortion. The effective superexchange 
between electrons on neighbouring rungs is different for 
these two clusters. Therefore we obtain an alternation of 
3 e ff along the ladder direction. We did not calculate the 
mlsffective superexchange for the V2i-V22-ladders, since it is 
strongly modified by charge ordering. |34| All results are 
found in table O . 

For the Vi-ladders we find a large alternation of J e // 
along the ladder with Sj « 0.25, larger by an order of 
magnitude than we had estimated from the change of the 
vanadium-oxygen bond lengths. The dimerisation of tz 
along the ladder is St L ~ 0.2 and therefore also much 
larger than what we estimated from the variation of the 
bond length above. We also note that Sj and St L are of 
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Table 11. Effective vanadium- vanadium hoppings t, on-site- 
energy e and superexchange J e // for the distorted phasegiven 
in eV. Denomination of hoppings as given in figure hi Su- 
perexchange 3 e ff has been calculated with U — 4 eV, and 
Vl = 0.344 as for the undistorted phase (see text). 



Hopping 


Layer 


a Layer b 


V2i-V22-ladders 


tfl,CO 


0.173 


0.171 


ti.CO 


0.052 


0.057 


to.co.i 


0.060 


0.061 


tu,CO,2 


0.084 


0.086 


tcV21 a ^21b 




0.015 


^ c 1 a V*2 2 6 




0.010 


tcV22a ^216 




0.005 


tcV22a ^22b 




0.004 


ey 2 i 


1.265 


1.250 


e V22 


1.421 


1.428 




Vi-ladders 


tit.SD 


0.179 


0.180 


t£,SD,l 


0.062 


0.068 




0.040 


0.046 


to,sr>,i 


0.070 


0.070 


tD,SD,2 


0.061 


0.061 


tc,SD,l 




0.007 


tc,SD,2 




0.006 


evi 


1.350 


1.338 




0.066 


0.072 


Je//,2 


0.038 


0.044 


<5j 


0.263 


0.247 




Interladder 


t/Ll 


-0.081 -0.072 


t/i2 


-0.130 -0.128 



similar size, whereas we would expect a factor of 2 differ- 
ence from the standard superexchange estimate. The lat- 
ter fact is, however, easily explained: We estimate the ef- 
fective superexchange between electrons on neighbouring 
rungs of the same ladder by looking at the singlet-triplet 
gap in such a two-rung cluster. This gap essentially de- 
scribes the kinetic energy gain associated with having the 
option of locating both electrons on a single site, if they 
are in a singlet state, whereas they cannot do so if they are 
in a triplet state. This means that all possible two-particle 
states contribute to the size of the singlet-triplet gap, in- 
cluding the states where both electrons are on a rung or 
sitting on the diagonal. Therefore the standard estimate 
using J = 4t 2 /U must fail and this is also the reason for 
our definition of J e //. This explains why Sj and 5t L are 
of similar size. 

The reason for the unexpectedly strong alternation 
of both J e ff and along the Vi-ladder of both lay- 
ers is more complicated. To understand this we take a 
closer look at how ti, is actually derived from the original 
vanadium-oxygen hopping terms and how these change 
under distortion. Consider hopping along the leg from 
vanadium site V to V via the intermediate Oi as shown 
in figure |L The matrix elements are given in table |l^. We 
observe that due to the signs of the hopping elements the 



Table 12. Vanadium-oxygen hopping along the leg given in 
eV. Nomenclature as in figure [| The first two lines give the 
hopping for the undistorted phase, the last four lines along the 
rungs of Via-ladder in the distorted phase. 



'Pxd X y ^Pyd X y ^Vzd X y 



undistorted phase 



VOi 


-0.8549 


0.6303 


-0.2470 


V'Oi 


0.8549 


0.6303 


0.2470 


distorted phase: 


layer a 




VOi 


-0.8486 


0.6758 


-0.2534 


V'Oi 


0.8486 


0.6758 


0.2534 


V0 3 


0.8676 


0.5860 


0.2358 


V"0 3 


-0.8676 


0.5860 


-0.2358 



effective hopping from V to V via the Oip y orbital coun- 
teracts the effective hopping via the Oip^ and the Oip z 
orbital. Thus, in a simple approximation, which neglects 
all other hoppings taking place in the system, the total 
effective hopping is 



t 2 t 2 t 2 

t vv , = E^l. + p y d *y p^v_ = _ .250 eV. 

e v — to £v — eo £y — eo 

_ (8) 
As we can also see from table [12], distortion affects the hop- 
pings via the different oxygen-p-orbitals in a significantly 
different way. For layer a the dimerisation of these hop- 
pings are approximately S x = 0.01, S y — 0.07, S z — 0.04. 
But whenever l^w J decreases, Lj „ increases and 

| tf\c xy | ' | r^y zcy | 

vice versa. Thus the total dimerisation is enhanced: With 
the same crude approximation as in equation (JsJ) we ob- 
tain tvv = -0.207 eV and tvv = -0.294 eV. This gives 
a dimerisation of the effective vanadium-vanadium hop- 
ping of St vv — 0.17. Thus we can explain the very large 
alternation of the superexchange J e // which we obtained 
for the Vi-ladders. Thus the main result of the section is 
that this alternation is almost completely due to the shift 
of the oxygen atoms on the leg of the Vi ladders in a and 
c direction, and is not caused by the small shift of the 
vanadium sites. 

We can also examine the behaviour of the hopping be- 
tween layers. For the spin-dimerized Vi ladders the size of 
the effective hopping t c ,SD changes slightly, but there ap- 
pears to be no significant alternation in c direction. This is 
different for the effective hopping between charge ordered 
ladders of neighbouring layers: Hopping is increased be- 
tween V21 sites and decreased between V22 sites. The size 
of the effective hopping t c ,co remains very small however, 
causing a superexchange coupling of less than 1 meV be- 
tween completely localized electrons on the V21 sites. We 
conclude that both in the undistorted and the distorted 
phase all important spin interactions take place within a 
vanadium-oxygen layer. 
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5 Summary 

In this article we report a structure determination by syn- 
chrotron X-ray diffraction of the low temperature phase of 
a'-NaV2C>5. As in [0,|l^] two slightly different vanadium 
layers are found to exist as well as two types of two-leg 
vanadium ladders within each layer. One of these ladders 
exhibits a zig-zag-type modulation and contains two in- 
equivalent vanadium sites V21 and V22- The other ladder 
contains only one vanadium site Vi . 

Based on the structure determination we perform a 
bond-valence calculation. Contrary to Q we find only two 
different valences: the V22 and the Vi sites experience 
an electron depletion towards a valence of v 4 5+l5 , while 
the charge on the V21 sites is increased towards a valence 
of V 4 - 5 ~ 315 . We compare our results with results of a 51 V- 
NMR measurement || and find them to be in qualitatively 
good agreement: Due to charge conservation an asymme- 
try both in the shift of the two new peaks from the original 
V 4,5 peak as well as an asymmetry in intensity are found. 

We then use a Slater-Koster treatment to find the hop- 
ping along the bonds and project out the V-0 and V- 
Na processes to obtain effective vanadium- vanadium hop- 
pings. The results compare well with LDA || as well as 
ab initio calculations [p9[ regarding relations between the 
different matrix elements on the ladder. When we cal- 
culate the effective superexchange between rungs for the 
undistorted phase using our matrix elements we find a 
value of 0.047 eV which compares well with experimen- 
tal results of 0.045 ... 0.048 eV. In the distorted phase 
we obtain a strong alternation of the superexchange with 
dimerisation Sj ss 0.25 in the Vi ladders along ladder di- 
rection. This dimerisation is nearly exclusively due to the 
shift of the oxygen atoms mediating the hopping along 
the legs of the ladders in a and c direction; the shift of 
the vanadium atoms is much smaller. Even though the 
Slater-Koster method is known to overestimate the effect 
of small distortions, the results show clearly, that a signif- 
icant alternation of the superexchange is taking place on 
the Vi-ladders. No alternation in b-direction is found for 
the V21-V22 ladders. No significant hopping or superex- 
change between layers is found both for the undistorted 
and the distorted phase. 

This leaves the question of the origin of the spin gap 
which has been observed in the material. According to 
our results for the valence of the different vanadium sites 
the picture of spin-singlet clusters proposed in |12| can- 
not be correct, since it relies heavily on having sites 
in the distorted phase. However, in 23 Na-NMR eight in- 
equivalent sodium sites have been detected [jl8|,[l0|. In our 
structure determination we obtain only six incquivalent 
sites. Four of these sites are located above and below the 
spin dimerised Vi -ladders, two of them are located above 
and below the charge-ordered V2i-V22-ladders. To have 
eight inequivalent sodium sites it would be sufficient to 
assume, that along the V2i-V22-ladders the shift of equiv- 
alent sodium sites in c direction alternates along the ladder 
direction. This would also immediately imply an alterna- 
tion of the superexchange on the charge-ordered ladders, 
since, as we have shown above, the hopping via the sodium 



site is crucial to explain the diagonal hopping between two 
rungs of the same ladder. Alternatively, such an alterna- 
tion of inequivalent sodium atoms can result from an al- 
ternating shift of the rung oxygens of the V2i-V22-ladders 
in b direction. This would also cause a superexchange al- 
ternation along the V2i-V22-ladders and one would also 
have eight inequivalent sodium sites in the system. Such 
shifts would explain in a simple manner both the observa- 
tion of 23 Na-NMR and the existence of a spin gap without 
disagreeing so far with the results of all X-ray structure de- 
terminations. Since the superexchange alternation seems 
to be much stronger on the Vi-ladders, this would imply 
a system with two significantly different spin gaps. 

Another explanation for the observation of eight in- 
equivalent sodium sites has been given in Jh}. The authors 
propose that a charge ordering takes place on all ladders. 
However, this implies having eight incquivalent sites pro- 
vided, one has a) an alternation of sodium sites along the 
ladder direction, thus implying superexchange alternation 
as above for the case of the V2i-V22-ladders, b) two in- 
equivalently charge ordered ladders. If one had the same 
amount of charge ordering on both ladders one would also 
expect the same sodium sites, since charge/ valence and 
bonds are strongly related. So, in such a case, four in- 
equivalent vanadium sites would have to be observed in 
51 V-NMR, which is not the case [§. 

In conclusion we have done a structure determina- 
tion of the low temperature phase of o/NaV2 05. We con- 
firm the results of [Q,[l2), having charge ordered V21-V22- 
ladders and disordered Vi-ladders. However, a bond va- 
lence calculation using our results finds only two instead 
of three significantly different valences: The valence of the 
Vi and the V22 sites turns out to be nearly the same. By 
making use of the Slater-Koster method we find a strong 
alternation of the superexchange coupling between rungs 
of the Vi-ladders. 

T. C. would like to thank Dr. H.-G. Krane for his help doing 
synchrotron X-ray measurements. 
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